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ON THE ZEROS OF RANDOM HARMONIC POLYNOMIALS: 
THE TRUNCATED MODEL 

ANTONIO LERARIO AND ERIK LUNDBERG 


Abstract. Motivated by Wilmshurst’s conjecture and more recent work of W. 
Li and A. Wei HU, we determine asymptotics for the number of zeros of random 
harmonic polynomials sampled from the truncated model, recently proposed 
by J. Hauenstein, D. Mehta, and the authors m- Our results confirm (and 
sharpen) their (3/2)—powerlaw conjecture [10] that had been formulated on the 
basis of computer experiments; this outcome is in contrast with that of the model 
studied in m- For the truncated model we also observe a phase-transition in 
the complex plane for the Kac-Rice density. 


1. Introduction 

A harmonic polynomial is a complex-valued harmonic function given by: 

(1) F{z) =p{z) + q{z), 

where p and q are polynomials of degree n and m (respectively). Let Np denote 
the number of zeros of F, that is, points z G C such that F{z) = 0. 

For n > m, we have the following bounds: 

n < Np < n^. 

The lower bound is based on the generalized argument principle and is sharp 
for each m and n. The upper bound follows from applying Bezout’s theorem to 
the real and imaginary parts of F{z) = 0 after noticing that the zeros are isolated, 
which was shown by Wilmshurst [29]. 

1.1. Wilmshurst’s conjecture. Wilmshurst made the conjecture that the Be- 
zout bound can be improved to a function that is linear in n for each hxed m, 
namely: 

(2) Np < 3n — 2 -|- m{m — 1) (Wilmshurst’s conjecture) 

This conjecture is stated in [23 Remark 2] (see also |2S] and |1]). 

For m = n — 1, the upper bound follows from Wilmshurst’s theorem [29], and 
examples were also given in [29] showing that this bound is sharp (shown indepen¬ 
dently in [2]). For m = 1, the upper bound was shown by Khavinson and Swiatek 
[IT] using anti-holomorphic dynamics. A proof of the Crofoot-Sarason conjecture 
given in [8] (cf. [3]) established that this bound is sharp. Counterexamples to 
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the case m = n — 3 were established analytically in [15] , and counterexamples for 
a broad range of (hnitely many) m and n were established in [10] using certihed 
numerics. On the other hand, we still expect, in the spirit of ([2]), that Np satishes 
an upper bound that is linear in n for m hxed; for instance, with S-Y. Lee, the 
authors conjectured in [151 Introduction] that Np < 2m{n — 1) + n. 


1.2. A probabilistic version of the problem. Given the high variability of the 
number of zeros Np^ it is natural to ask the following. 


Question 1. What is the expectation EiNp of the number of zeros of a random 
harmonic polynomial? 


This question was asked and answered by W. Li and A. Wei in HZI, in the case 
when p and q are independently sampled from the complex Kostlan ensemble: 

n m 

(3) p(z) = ^afcZ^ q{z) = '^bkZ^, 

k=0 k=0 


where Ok and bk are independent centered complex Gaussians with KajOk = djk{'^) 
and Ebjbk = 

The choices of p and g in ([3]) lead to the following asymptotics (as n —)■ oo): 

(4) ~ I W"' ,„ 

1^ n, when m = an + o[n) with U < a < 1, 


Notice that when m = an the average number of zeros is asymptotically the 
fewest possible. This seems to suggest that, on average, an even stronger form 
of Wilmshurst’s conjecture ([2]) holds. However, caution is needed here, and the 
dichotomy in (|1|) dissolves after choosing a dehnition of “random” in which the 
coefficients of p and q are more comparable in modulus (see Theorem [1] below). 

In the model ([3]), where the coefficients of p are asymptotically much larger 
in modulus than q when m = an, F tends to resemble an analytic polynomial 
and asymptotically obeys the fundamental theorem of algebra. In order to make q 
more comparable to p, an alternative model (referred to as the “truncated model”) 
was proposed in m where the variances were replaced by in the dehnition 
fl51) of q, while still choosing m as the upper limit in the summation (see dehnition 
fl51) below). For the truncated model, computer experiments performed in [TU] led 
to a conjecture that the expectation ENp has a (3/2)—powerlaw growth whenever 
m = an for all 0 < a < 1. Here, we prove (and sharpen) this conjecture, see 
Theorem [1] below. 

Note that we do not consider here the case m = 0 of random complex analytic 
polynomials, where we would have Np = n almost surely (by the fundamental 
theorem of algebra). Yet, it is still interesting in that case to study the location of 
zeros; we refer the reader to Edelman and Kostlan’s paper [?[ Sec. 8] and to the 
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recent work of Zeitouni and Zelditch [30] establishing a large deviation principle 
for the location of the zeros of a random analytic polynomial. 


1.3. Asymptotics for the truncated model. We revisit Question [T] while sam¬ 
pling F{z) = p{z) + q{z) randomly from the truncated model, i.e., 

n m 

(5) p{z) = ^akZ^, q{z) = ^bkZ^, 

k=0 k=0 

where ak and bk are independent centered complex Gaussians with EojOfc = 6jk(j) 
and Ebjbk = bjkQ). 


Theorem 1. Let F{z) = Pn{z) + qm{z) be a random polynomial from the truncated 
model. For m = an with 0 < a < 1, the expectation KNp of the number of zeros 
of F[z) satisfies the following asymptotic (as n ^ oo) 


where Ca is given by 


( 6 ) 



On the other hand, when n —)■ oo with m fixed we have ENp ~ n. 


Our methods can be used to describe asymptotics for the Kac-Rice density 
(providing the expected number of zeros over a prescribed region). We notice a 
phase-transition in this pointwise asymptotic, and the leading contribution 
is completely accounted for by zeros that are located within a critical distance 
from the origin, see Section 13.31 

Note that as a —)■ 1, Cq, 7r/4, in agreement with [T^ Thm. 1.1]. 

An interesting aspect of harmonic polynomials is that, unlike analytic poly¬ 
nomials, the function F{z) = p{z) + q{z) can reverse orientation. The orienta¬ 
tion of F can be determined by the sign of the Jacobian determinant Jp{z) = 
\p'{z)\^ — \q'{z)\^. Let A_|_ denote the number of zeros for which F is orientation¬ 
preserving (i.e., Jp < 0) and A_ denote the number of zeros where F is orientation- 
reversing {.Jp > 0). 

Using a standard application of the generalized argument principle, we then 
notice the following corollary of Theorem [H showing that orientation-reversing 
zeros are asymptotically as common as orientation-preserving ones. 

Corollary 2. For m = an with 0 < a < 1, we have EA+ ~ EA_ ~ . 

Proof. Almost surely we have Np = Np + A_ (the presence of singular zeros 
is a probability zero event). By topological degree theory (or the generalized 
argument principle [5]) the difference Np — N_ is given by the winding number 
of F along a sufficiently large circle. Moreover, since the z” term dominates, the 
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Figure 1. A portion of a random critical lemniscate (the critical 
set of a random harmonic polynomial) with m = n = 100. Plotted 
in the region {z G C : |3fJz| < 1, IQ'zl < 1}. For m = n the truncated 
model coincides with the Li-Wei model. 

winding number is n, and so we have N+ = N- + n. Theorem [T] then implies that 
EN+ ~ EN_ ~ f n3/2. □ 

The coexistence of many zeros of opposite orientation suggests that the Jacobian 
of F changes sign wildly throughout the complex plane (or otherwise that there is 
a high level of “condensation” of zeros into regions of common orientation). Taking 
this point into consideration, we conclude the introduction by posing the problem 
of investigating the topology of the orientation-reversing set := { 2 ; G C : 
|p'(z)| < |g'( 2 ;)|}. It follows from applying the maximum principle to the harmonic 
function log \p\z) \ — log \q\z)\ that each connected component of contains at 
least one critical point of p. This implies that has at most n — 1 connected 
components. What can be said about the average number of components of f2_? 
The critical set (the boundary of f2_) is depicted in £gure[T]for a random sample 
with m = n = 100. Note that the critical set is a random rational lemniscate, 



(7) 


similar to the random lemniscates studied recently by the authors [16]; the only 
difference is that in the model studied in [161 Sec. 1.2], the numerator p and 
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denominator q of the rational function appear without differentiation. Based on 
the results in [TB] we conjecture that when m = n ^ oo the average number of 
connected components of the random critical lemniscate ([7]) grows linearly (the 
maximum rate possible). 

1.4. Outline. In Section O we provide an exact formula for the average number 
of zeros for the truncated model. This is derived from a slight modihcation of 
HZI. The asymptotics stated in Theorem [1] are proved in Section |3l The proof 
uses the dominated convergence theorem after factoring out Establishing a 
dominating function requires several elementary estimates, and determining the 
pointwise limit of the integrand requires asymptotics for a truncated binomial sum. 
Such asymptotics are provided in Lemma IU and the proof of Lemma H] is given in 
the separate Section 01 The proof uses both forms of Laplace’s asymptotic method 
m Sec. 3.3, 3.4]: namely the case of an interior maximum (saddle-point) as well 
as the case of an end-point maximum. The presence of both cases is responsible 
for the phase-transition in the Kac-Rice density mentioned above (see also Section 

E31). 

Acknowledgement. We wish to thank Seung-Yeop Lee for inspiring discussions 
and insightful suggestions and the anonymous referees for their helpful remarks. 


2. An exact formula for ENf 

Let Pm,n{x) := denote the binomial expansion of (1-1-x)” truncated 

at degree m. 


Theorem 3. The expectation ENpiT) of the number of zeros of Fn,m{,z) = Pn{z) + 
qm{z) on a domain T G C is given by: 


( 8 ) 


ENf{T) 


1 f 1 RI + RI-2RI, 
7rJT\z\^Ry{R, + R,y-4Rl, ^ 


where dA{z) denotes the Lebesgue measure on the plane, and 

R,2 = n^\z\\l + \z\^r-^P^.,,n-i{\z\^), 

Rs = (1 + kl^)” + Pm,n{\z\^), 

Ri = R:i{n^\z\* + n\z\^){l + + \z??"~\ 

R 2 = R3[n‘^\z\^Pm-2,n-2{\z\‘^) + n\z\‘^ Pm-l,n-2{\z\^)] - n^\z\‘^[Pm-l,n-l{\z\‘^)? ■ 


Note: The analogous statements contained in [171 Thm. 1.1, Thm. 4.1] contain 
a little ambiguity. In fact the authors use the Kac-Rice formula for the harmonic 
function already in polar coordinates, thus viewing it as a random held dehned 
over [0, 27r) x (0, 00 ) and with values in R^. In particular [TTl Equation (1.1)] 
should either be modihed with \z\^ instead of |z| (and da{z) is still the Lebesgue 
measure on the complex plane) or the integration should be performed over the 
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image of T C C under the polar change of coordinates (and in this case \z\ = p). 
In other words, denoting by -0 : C!\{0} —)■ (0, 27r) x (0, cxd) the polar change of 
coordinates, the right expression for m Equation (1.1)] is: 


(9) r 


+ r| - 2 r 


12 


TT Jt ^3\/(^i + '1^2^ - 4ri2 


:da{z) = — 


rj + rl - 2r 


12 


TT Jp(T) rl \/(r7Tr^P^^4^ 


-pdpdO 


This ambiguity is no longer present in their asymptotic analysis. 


Proof. We follow closely the lines of the proof given in [T^ Thm. 1.1, Thm. 4.1], 
adjusting certain computations as needed. Also, we simplify the hrst part of their 
proof by not switching to polar coordinates while obtaining equations ffTTll and 
flTSl) below. 

Applying the Kac-Rice formula (restated in [T71 Lemma 2 . 1 ]), we have: 


(10) ENp{T) = jE{\detJF{z)\\F{z)=0)p{0;z)dA{z), 

where, for each z, p(s; z) is the probability density function of the random variable 
s = F{z). 

The modulus of the Jacobian determinant of F{z) = p{z) + q{z) is given by (see 
|ni Sec. 1.2]) 

(11) |*(j)| = ||pn)|= - W(z)f\ = F\\zp'(z)r- - \zq'(zp \, 
and hence we have 

(12) E (I det Jf{z)\\F{z) = O) = [\ul — ul + vl — = OWa = O) i 

where, for j = 1, 2, 3, the expressions Uj, Vj are given by 


Ml = 3? kttkZ^, 

k=0 

m 

k=0 


Ml = A kokZ^, 

k=0 

m 

V2 = ^y^^kbkz’", 

k=0 


U3 


[Pniz) + qmiz) 



Then letting (f/i, U 2 , Vi, V 2 ) denote the Gaussian vector that has the distribution 
of (mi, M 25 '^^ 1 ^ 2 ) under the (linear) condition M 3 = 0 ,M 3 = 0 , we have 

(13) E {\ul -ul + vf - vl\\u3 = 0,M3 = 0) = E {\U^ - 17| + _ vi\) . 

The covariance matrix R of (f/i, U 2 , R, R) is given by [26l p. 30] 

(14) R = C - 
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where 


^2x2 = COv(M3,t;3), 

^4x2 = COv((mi, M2, Vi, V 2 ), (^ 3 , V 3 )), 
C4x4 = COv(mi,M2,Pi,P2)- 


First we compute 


Eul = EU 32 = -J2 


k=0 


\2k 


+1E f = ^(1 + '•■I')” + 


fc =0 


Euiug = Euir;3 = \ \z\^^ = \ 


k=0 


k=0 


E«; = E«J = ;; E U 


1 m / \ 

EM 2 M 3 = EU 2 P 3 = 9 l) 1 ^ 1 ^^ = 9 ’^kP^m-l,n-l(kn, 


k=0 


Eul = Ep2 = = ^(,n‘^\z\^ + n\z\‘^)Pm-2,n-2{\z\‘^), 


k=0 


and hence 


^2x2 — 


(i + izn-+p^,„(izn A 0 


/ n\zWl + \z\^)^-^ 


Bax2 — 


0 ir 

0 \ 

n\z\‘^Pm-l,n-l{\z\‘^) 0 

0 n|z|2(l +|z|2)-i 

y 0 n\z\‘^Pm-l,n-l{\z\‘^) J 


Caxa = -diag((A|2;p + n)(l + \z\'^)'^ ^ {n^\z\^ + n)P„,_2,n-2{\z\‘^), 

(A|z |2 + n)(l + |z| 2 )"- 2 , + n)P^_2,.-2(kP)). 


From these we compute 


R 


4x4 — 


1 


/ 


—-R12 

0 

0 

\ 


-Ri2 

R2 

0 

0 



0 

0 

Ri 

—Ri2 

V 

0 

0 

— -R12 

R2 

/ 
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where 

R,2 = nM\^ + \z?T~"Pm-i,n-i{\z?). 

/?3 = (i+krr+^™,n(kr), 

Ri = R‘i{n^\z\^ + n\z\^){l + \z\^Y~'^ - n^\z\^{l + \z 
R 2 = R3[n‘^\z\*Pm-2,n-2{\z\‘^) + Ti 12;| ^Pm_i,n-2 ( ^ H] “ \z\'^[Pm-l,n-l{\z\‘^)f ■ 

Applying [TTl Cor. 2.1], we obtain: 


|2\2n-2 


(15) 


E\u^ -ui + v^ -vi\ = 


1 Rl + Rl- 2P2 


12 


R3^{Ri + R2)^-ARl^' 

For each hxed z, the complex Gaussian s = F[z) has probability density func¬ 


tion 


and in particular 


p{s-z) = —-exp{-|s|VP 3 }, 
ttPs 


p(0;^) = 


1 


ttP.s 


Applying this along with equations (IT^ . (IT^ . and (1T5|1 to the Kac-Rice formula 
m we obtain the desired result ([ 8 ]). □ 

3. Proof of Theorem [U 

3.1. The case when m = an. Applying Theorem |3] with Np := Np{C), switch¬ 
ing to polar coordinates r = \z\, dA{z) = rdrdO, and integrating out the angular 
variable 6, we are left with: 

ENp = 2 

where 


f°° 1 -|- 02 — 2a^2 

'0 a|\/(ai 02 )^ - 40^2 


-.dr, 


012 = nV(l -h ^Pm-i,n-i{r^), 

03 = (1 + + Pm,n{Y), 

Oi = -I- nr^)(l -|- — n^r‘^(l -|- 

02 = 03[nV‘Pm_2,n-2(r^) nr^Pr„-i,n-2(r^)] - n^rYPm.-i,n-i{Y)f . 
Factoring (l-|-r^)‘^’^“^ from the numerator and (l-|-r^)‘^’^“^ from the denominator, 
we have: 

1 + 


(16) 


ENp = 2n^R 


Iq nR^r{l + rYHlYibi + b2Y-^bl 

^Rm—ljn—lY ) 


-.dr, 


12 


612 = nr 


(1 -F r2)"-i ’ 


where 
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= 1 + 


m,n\' ) 


bo = ho 


nr 


4 -^ 771 — 2 , 71 — 2 ( 1 " 


{l + r'^Y' 
hi = h^[nr^ + r^] — nr^, 

2-^ 771 —l,n—2 ) 


+ r 


nr 


Pm—l,n—l{j' ) 

_ (1 + r 2 )’"-i _ 


-f- 2 _j_ ^2^n—2 

We will apply Lebesgue’s dominated convergence theorem to take the limit of 
the integral appearing in (1T6|) . The following claim implies that the sequence of 
integrands in fITB]) is bounded by a single integrable function. 

Claim: 

bl + hl- 26^2 r- 3^ 

- , = 0 (vnr ), as n -)■ 00 . 

hl^WpW^2 

Proof of Claim. First we note that 0102 > ai 2 - This is by the Cauchy-Schwarz 
inequality, since it follows from the proof of Theorem [3] that oi = Et/f, 02 = Et/|, 
and ai 2 = WJ 1 U 2 , where Ui and U 2 are Gaussian random variables. 

This implies that 6162 > ^ 12 - Since 63 > 1, we have: 

hl + hl- 26^2 


bW{hi + h2y-Pbl^ 


<\ hl + hl- 2bl 


12 


= ^J{hl-h2f + 2{hlh2-hl^) 

^ \/{hi — 62)^ + \J 2(& i &2 ~ &12) 
= 161 — 621 + &1&2 


G2 

Thus, it suffices to show that 

(17) 61 - 62 = 0(\/u + r^), 
and 

(18) 6162 — ^12 = 0 (ur®). 

Let gm,n := Then, 63 = 1 + and we have: 

61 = (1 + qm,n){nr'^ + r^) - nr^, 

and 

62 = (1 + qm,n) {nr^qm- 2 ,n -2 + r‘^qm-l,n- 2 ) - „_ 3 . 

These lead to: 

bi-b 2 = nr^ ((1 + qm,n) 0 - “ ^771-2,71-2) - (1 - g^_i,„_i)) 
+ (1 + qm,n)^‘^ (1 ~ qm-l,n- 2 ) ■ 
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The term (l + gm,n)i’^ (1 — <?m-i,n- 2 ) is bounded by 2r^. We consider the remaining 
term nr‘^ ((1 + gm,n)(l - <?m- 2 ,n- 2 ) - (1 - which can be rewritten as: 

flT ((1 T (lm,n){,Q.m—l,n—l Qm—2,n—2^ T (1 (lm—l,n—l)iSlm,n l,n—l)) 

^TIT (lm—2,n—2} T i,Qm,n Qm—l,n—l)') 


=nr 


<nr^ 


<3n 


n — 2 


^2{m-l) 


m — 1/ (1 + 


+ 


77, — 1 


^2m 


n — 2 


+ 


m — 1 

fl — \\ j,2{m+l) 


n — 1 
m 


m / (1 + 

^2(m-l) 


(1 + 


m y (1 + ’ 

where we have used the identity 


(19) 


Q.m,n Qm—l,n—l 


n — 1 


^2m 


m / (1 + 


which can be seen as follows 

Er.oG)’-“-(i+OEr.-o 




(1 + 

Er.ih(;)-(T), a-a -■ 


n 


(1 + r'^y 

r^2m 


2m 


m / (1 + r‘^y 


m / (1 + r^)"- 

Applying the hrst derivative test to over the interval x > 0 we hnd that 

the maximum occurs at x = Thus, we have: 

n—m—2 ’ 

^2(m+l) 

_ ‘1^ I I _2: 

m J ( 

/ ^n—m—2^ 

- 1\ (m + - m - 2 )^-^-^ 


3n 


n — 1 


m / (1 + 


< 3n 


< 3n 


m 


{n - I)*"-! 


= 3n- 


m + 1 f n — l\{m + — m — 2) 


< Cn 


n — 1 — m \m + 1 
m + 1 


n—m—2 


n 


n — 1 — m]J {m + l){n — m — 2) 


{n - 1)^-1 

= 0(n'/% 


where we have used Stirling’s approximation while recalling that m = an. 
This establishes flT7|) . 

Next we consider 6162 — &i 2 - 
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We have: 

blb 2 - [(1 + qm,n)qm- 2 ,n -2 “ 5^-1,n-l] “ 


I /I I \ 2 f I Pm—l,n—2{j' ) \ 

+ (1 + qm,n)r ^2^n-2 j ' 


Part of this can be estimated as follows: 

(1 + q^,n)r‘^ {b 2 + < 2r^(2(nr^ + r^) + 1) = 0(nr®). 

Since 6162 — &12 > 0, in order to prove flTSll . it is enongh to show that for the 
remaining terms we have: 

1’ {qm,n [(1 + qm,n)qm-2,n—2 ~ ~ — 0- 

We notice that 

{qm,n [(1 + qm,n)qm-2,n-2 ~ ~ 

(Q'm,n(l “1“ qm,n)qm—2,n—2 (1 “1“ ?m,n)l,n—l) 

= (1 + qm,n) {qm,nqm—2,n—2 l) • 

We will show that qm,nqm- 2 ,n -2 - < O' 

Using again the identity (IT^ . we have: 

'n-2\ 


qm,nqm—2,n—2 l,n—1 qm,n 1 l,n—1 


m-iy (l + r2)’^-iy 
n-2\ 
m - 

n-2\ 
m — 1/ (1 + 


9m—l,n—1 


m 


— 9m—l,n—1 I9m,n 9m—l.n—l) 9m,n I -i I /-i , o\^ i 

\m — 1 / (1 + 


n — 

j (1 + r2)" ■ 


r ,n -1 


zr ^ -Pm-l,n-l(L ) -Pm,n(L ) 


(1 + r2)2”--l [ 77i 


Finally, we have: 

oU — 1 


m 


- Pm,n(L^) = “1 + X] 


i=i 


n — 1 /n — 1 





-2i 


m \j — 

and we see that each coefficient ~ (j) = (jZi) “ f) negative. 


□ 


Having jnstihed an application of Lebesgne’s dominated convergence theorem, 
we hnd the pointwise limit of the integrand in (fT 6 |) nsing the following asymptotic 
(whose proof is given in Section 0]) . 
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Lemma 4. Let x > 0. For all 0 < a < 1, we have (as n ^ oo with m = an): 

Pm,n{x) _ [ l + 0(l/n), 0<X<j^, 

(1 + x)” “ \ 0(exp{-cn}), x > 


According to this asymptotic, for the integrand in Equation ffT 6 |) 

converges to zero, and for 0 < we see that b 2 = bi{l + 0{l/n)), and 

nF%l^{bi + b2y -Abl^ nP^bl ^ ^ ^ 2 ' 


Thus, we have 


Nf ~ 


’\/“/(l-“) ^2 


where 


(1 + r2)2 
-y/a/il-a) ^2 


dr = n^Pca, 


Cry - 


rdr. 


(1 + r2)2 

In order to determine Cq, we make the change of variable r = tan( 6 '),(ir = 
sec‘^{9)d9: 

l‘^/a/(l-a) ^2 rA 

-dr = / sin^( 6 ')d 6 ', A = arctan 


(1 + r2)2 


a 


\ — a 


Thus, we have 


Ca = ^ ^arctan (y\j ^ ^ j “ \/oi{l - a)^ . 

This completes the proof of Theorem [1] in the case that m = an with 0 < a < 1. 


3.2. The case when n —)■ oo with m fixed. This case is simpler and does not 
require Lemma 01 Omitting the details, we hnd that & 2 ,&i 2 ; converge to zero, 63 
converges to 1 , and 


Np ~ 



r 

(1 + r2)2 


dr = n. 


3.3. Asymptotics of the Kac-Rice density. Consider again the case when 
m = cm, with 0 < a < 1. Above, we have factored out n^P from the Kac-Rice 
density in order to apply the dominated convergence theorem, but Lemma 0] can 
also be used to hnd the pointwise asymptotic. The Kac-Rice density is asymptotic 
(as n -A- 00 ) to 27 rq-|-|i^ 2)2 for \z\ < ■\/a/(l — a), and it is asymptotic to 
for \z\ > ^Ja/{l — a). Thus, the leading contribution of zeros are located within 
the distance -^ 0/(1 — a) from the origin. This critical radius originates in the 
proof of Lemma 0] (based on Laplace’s method), see Cases 1 and 2 in Section 01 
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4. Proof of Lemma 0] using Laplace’s method 
The following formula is provided in [201 Lemma 1]. 

Lemma 5. For 0 < m < n — 1 


( 20 ) 


(I + X)^ 


n 


m 


{n — m) 


F/{x+i) 


u^{l-u) 


n—m—1 


du. 


We apply Laplace’s method to derive Lemma 0] from Lemma [5l Rewriting the 
integrand, we have for m = an: 

J x/{x+l) 

where h{u) = [a log(M) + (1 — a) log(l — m)], and g{u) = (1 — u)~^. 

Case 1: When x/{x + 1) < a, h{u) achieves its maximum at m = a, the unique 
solution of the saddle-point equation: 

h' {u) = a/u — (1 — «)/(! — m) = 0. 

Applying Laplace’s method [ISl Sec. 3.4], we have: 


I x/{x+l) 




27r 


(l + 0(n-^)) 


-nh"{a) 

Applying Stirling’s approximation, we have: 


( 21 ) 


in — m) = 


n 


{l + 0{n ^)) . 


— a)d ^-^27ra(l — a) 

Combining these results into (l20jl . we hnd: 

LiAA = (1 + o(n-')), 

Case 2: When x/{x -|- 1) > a, the saddle-point m = a is outside of the interval 
of integration, and h{u) instead achieves its maximum at the left end-point u = 
x/{x + 1). We thus have (by the alternative form of Laplace’s method [191 Sec. 
3.3]): 

Jx/ix+i) -nh'{x/{x + l))^ ^ 


X 


x+1 


an+l 


X -I- 1 


(1—Q;)n—1 


n(x(l — a) — a) 
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Combining this with (ET]), we see that 



(1 + a;)” 


r\j 


ci{x, a) 


^—C2{x,a)n 



5. Concluding remarks 

We have shown that the average number of zeros of a random harmonic polyno¬ 
mial sampled from the truncated model has order when m is a hxed fraction 
of n and grows linearly in n when m is hxed. In comparison with the Li-Wei model 
m Thm. 1.1], this behavior seems more indicative of (a probabilistic version) of 
Wilmshurst’s conjecture. 

Extending the above-mentioned breakthrough [T3], Khavinson and Neumann 
[12] used anti-holomorphic dynamics to count zeros of rational harmonic func¬ 
tions of the form r{z) + z, giving a complete solution to astronomer S-H. Rhie’s 
conjecture [2l| in gravitational lensing. For further discussion and related results, 
see usiniE]. In order to model stochastic gravitational lensing, the zeros of 
random harmonic functions were studied by A. Wei in his thesis [271 Ch. 3] and 
by Fetters, Rider, and Teguia [221 E3] • 

Is the variance of Np asymptotically proportional to the mean? Computer 
experiments in [10] suggest that the answer is (perhaps surprisingly, cf. [9]) “no”, 
and that the variance instead has order n^. 
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